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Abstract 

We investigate the tunnelling of spin 1/2 particles through event hori- 
zons. We first apply the tunnelling method to Rindler spacetime and 
obtain the Unruh temperature. We then apply fermion tunnelling to a 
general non-rotating black hole metric and show that the Hawking tem- 
perature is recovered. 

1 Introduction 

In recent years, a semi-classical method of modeling Hawking radiation as a 
tunneling effect has been developed and has garnered a lot of interest [T]-[20]. 
The earliest work with black hole tunnelling was done by Kraus and Wilczek 
[1] , an approach that was subsequently refined by various researchers [2 [3l [4] . 
From this emerged an alternative way of understanding black hole radiation. In 
particular one can calculate the Hawking temperature in a manner independent 
of traditional Wick Rotation methods or Hawking's original method of modelling 
gravitational collapse [21] . Tunnelling provides not only a useful verification of 
thermodynamic properties of black holes but also an alternate conceptual means 
for understanding the underlying physical process of black hole radiation. It 
has been shown to be very robust, having been successfully applied to a wide 
variety of exotic spacetimes such as Kerr and Kerr-Newmann cases [SI El [12] , 
black rings [1^ , the 3-dimensional BTZ black hole [3 [11] , Vaidya [16] , other 
dynamical black holes [TTj, Taub-NUT spacetimes [T^, and Godel spacetimes 
[20] . Tunnelling methods have even been applied to horizons that are not black 
hole horizons, such as Rindler Spacetimes [1|,[I1] and it has been shown the 
Unruh temperature [12] is in fact recovered. 

In general the tunnelling methods involve calculating the imaginary part of 
the action for the (classically forbidden) process of s-wave emission across the 
horizon, which in turn is related to the Boltzmann factor for emission at the 
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Hawking temperature. There are two different approaches that are used to 
calculate the imaginary part of the action for the emitted particle. The first 
black hole tunnelling method developed was the Null Geodesic Method used 
by Parikh and Wilczek U which followed from the work of Kraus and Wilczek 
[T] . The other approach to black hole tunnelling is the Hamilton- Jacobi Ansatz 
used by Agheben et al which is an extension of the complex path analysis of 
Padmanabhan et al [4]. Both of these approaches to tunnelling use the fact 
that the WKB approximation of the tunneling probability for the classically 
forbidden trajectory from inside to outside the horizon is given by: 

r oc exp(-2Im/) (1) 

where / is the classical action of the trajectory to leading order in h (here set 
equal to unity). Where these two methods differ is in how the action is calcu- 
lated. For the Null Geodesic method the only part of the action that contributes 
an imaginary term is J^°'^* prdr, where pr'is the momentum of the emitted null 
s-wave. Then by using Hamilton's equation and knowledge of the null geodesies 
it is possible to calculate the imaginary part of the action. For the Hamilton- 
Jacobi ansatz it is assumed that the action of the emitted scalar particle satisfies 
the relativistic Hamilton- Jacobi equation. From the symmetries of the metric 
one picks an appropriate ansatz for the form of the action and plugs it into the 
Relativistic Hamilton-Jacobi Equation to solve. (For a detailed comparison of 
the Hamilton-Jacobi Ansatz and Null-Geodesic methods see [H]). 

Since a black hole has a well defined temperature it should radiate all types 
of particles like a black body at that temperature (ignoring grey body effects). 
The emission spectrum therefore is expected to contain particles of all spins; 
the implications of this expectation were studied 30 years ago [23] . However 
application of tunnelling methods themselves to date have only involved scalar 
particles. Specifically there is no other black hole tunnelling calculation (to the 
best of our knowledge) that models fermions tunnelling from the black hole. 
In fact comparatively little has been done for fermion radiation for black holes. 
The Hawking temperature for fermion radiation has been calculated for 2d black 
holes [21] using the Bogoliubov transformation and more recently was calculated 
for evaporating black holes using a technique called the generalized tortoise 
coordinate transformation (GTCT) [lS]-[57]. The latter result is interesting 
because there is a contribution to the fermion emission probability due to a 
coupling effect between the spin of the emitted fermion and the acceleration of 
the Kinnersley black hole. From this one may infer that when fermions are 
emitted from rotating black holes that will be a coupling between the spin of 
the fermion and angular momentum of the rotating black hole present in the 
tunnelling probability. 

In this paper we extend the tunnelling method to model spin 1/2 particle 
emission from non-rotating black holes. In order to do this we will follow an 
analogous approach to the original approach used by Padmanabhan et al [4]. 
The Hamilton Jacobi ansatz emerged from an application of the WKB approx- 
imation to the Klein Gordon equation. We will start by reviewing this general 
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calculation, and then apply a WKB approximation to the Dirac Equation. We 
consider Rindler spacetime first and confirm that the Unruh temperature is re- 
covered. Insofar as fermionic vacua are distinct from bosonic vacua and can 
lead to distinct physical results [2H|, this result is non-trivial. We then extend 
this technique to general 4-D black hole metric and show the Hawking temper- 
ature is recovered. We illustrate this result in several coordinate systems - 
Schwarzschild, Painleve, and Kruskal - to demonstrate that the result is inde- 
pendent of this choice. This last system is particularly interesting since it has no 
coordinate singularities at the horizon. That we obtain the expected Hawking 
temperature indicates that tunnelling can be understood as a bona-fide physical 
phenomenon. 

One of the assumptions of our semi-classical calculation is to neglect any 
change of angular momentum of the black hole due to the spin of the emitted 
particle. For zero-angular momentum black holes with mass much larger than 
the Planck mass this is a good approximation. Furthermore, statistically parti- 
cles of opposite spin will be emitted in equal numbers, yielding no net change 
in the angular momentum of the black hole (although second-order statistical 
fluctuations will be present). We confirm that spin 1/2 fermions are also emitted 
at the Hawking Temperature. This final result, while not surprising, furnishes 
an important confirmation of the robustness of the tunnelling approach. 

2 Review of the Hamilton-Jacobi Ansatz 

We will consider a general (non-extremal) black hole metric of the form: 

ds^ = ~f(r)dt^ + —- + C(r)h,,dx'dx^ (2) 
9{r) 

The Klein Gordon equation for a scalar field is: 

9 

777 

Applying the WKB approximation by assuming an ansatz of the form 
(t){t, r, x') = exp[^I{t, r,x')+ h (t, r,x') + 0{h)] 

and then inserting this back into the Klein Gordon equation we get the usual 
result of the Hamilton-Jacobi equation to the lowest order in h: 

- [g^'^d^Id,! + m^] + 0{h) = 

(obtained after dividing by the exponential term and multiplying by h^). 
For our metric the Hamilton-Jacobi equation is explicitly 
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for the black hole metric ^ where we neglect the effects of the self-gravitation 
of the particle. There exists a solution of the form 

I = -Et + W{r) + J{x') + K (4) 

where 

dtI = ~E, drl^W'ir), dj^,h 

and K and the J^'s are constant {K can be complex). Since dt is the timelike 
killing vector for this coordinate system, E is the energy. Solving for W{r) yields 



W^{r) = ± / . ^ Je^ f{r){m^ + (5) 

since the equation was quadratic in terms of W{r). One solution corresponds to 
scalar particles moving away from the black hole (i.e. + outgoing) and the other 
solution corresponds to particles moving toward the black hole (i.e. - incoming). 
Imaginary parts of the action can only come due the pole at the horizon or from 
the imaginary part of K. The probabilities of crossing the horizon each way are 
proportional to 

2 2 
Prob[out] cx exp[— - Im/] = exp[— -(Im + Imif)] (6) 

2 2 

Probfm] cx expf Im/]=exp[ (ImVF- + ImX)] (7) 

h n 

To ensure that the probability is normalized so that any incoming parti- 
cles crossing the horizon have a 100% chance of entering the black hole we set 
ImK = — ImT4^_ and since W+ — —W- this implies that the probability of a 
particle tunnelling from inside to outside the horizon is: 

4 

r cx exp[--ImTy+] (8) 

Henceforth we set H to unity and also drop the subscript from W. 
Integrating around the pole at the horizon leads to the result [H] 

W= , (9) 

where the imaginary part of W is now manifest. This leads to a tunnelling 
probability of: 

47r 

r = cxp[ E\ 

and implies the usual Hawking temperature of: 



^ (10) 

It can be shown T9] that the proper Hawking temperature is recovered for 
multiple choices of the form of the metric for the same black hole. 
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3 Spin 1/2 particles and Rindler Space 



We first consider the Rindler spacetime, for which the tunnelHng calculation of 
a scalar field has shown that the Unruh temperature [H] is recovered. 

We will only show the calculation explicitly for spin up case; the final result 
is also the same for the spin down case as can be easily shown using the methods 
described below. Due to the statistical nature of the heat bath we assume that 
no angular momentum is imparted to the accelerating detector (i.e. on average 
there are as many spin up particles as spin down particles detected). The 
fermionic heat bath as seen by accelerated observers has many applications, 
such as understanding the effects of acceleration on entanglement [?] . 

We will use the following metric for Rindler spacetime 



9(z) 

f{z) = a^z^-\ 
a^z'^ - 1 



so chosen for its convenience in extending the technique to normal black holes. 



The Dirac equation is: 



where: 



vn 

i-i^'D^i, + -^V- = (11) 



= -AT. A 



4 



The 7'^ matrices satisfy {7'', 7^^} = 2f/'^'^ x 1. There are many different ways 
to choose the 7'' matrices and we will use the following chiral form: 



7* 



1/01 

-1 









a2 
^2 



0-3 

CT^ 
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where the a' s are simply the Pauh Sigma Matrices: 



1 

1 

~i 

1 

1 

-1 



and are the eigenvectors of . Note that 



-1 
1 



is the resulting 7^ matrix. 

Measuring spin in the z-direction (i.e. the direction of the accelerating ob- 
server) we employ the following ansatz for the Dirac field, respectively corre- 
sponding to the spin up and spin down cases: 



A(t,x,y,z)^^ 
B{t,x,y,z)^^ 

A{t,x,y,z) 


B{t,x,y, z) 


C{t,x,y,z)(,y 
D{t,x,y,z)^^ 



Cit,x,y,z) 


D{t,x,y, z) 



exp 



exp 



exp 



exp 



^I'\{t,x,y,z) 



:Il{t,x,y,z) 



^Iiit,x,y,z) 



:Il{t,x,y,z) 



(12) 



(13) 



In order to apply the WKB approximation we insert the ansatz (|12p for spin 
up particles into the Dirac Equation. Dividing by the exponential term and 
multiplying by h the resulting equations to leading order in h are 



B 



1 



dtl^ + yM^dJ^j +Am 

-B{d^I^ +idyl^) 



-A {dxl\ + idyl^ 



(14) 
(15) 
(16) 
(17) 
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Note that although A, B are not constant, their derivatives ~ and the compo- 
nents Vl^ - are all of order 0{li) and so can be neglected to lowest order in 
WKB. 

When m 7^ equations (|14p and (|16p couple whereas when m — Q they 
decouple. We employ the ansatz 

I^ = -Et + W{z) + P{x,y) (18) 

and insert it into equations (|14m7p 

-E 



B^^-j= + ^)W'{z)ymA = (19) 
-B{P,+iPy) = (20) 



E 



-A^^-^= + ^)W'{z)ymB = (21) 
-A{P.,+iPy) - (22) 

where we consider only the positive frequency contributions without loss of 
generality. Equations (PO)) and both yield {Px + iPy) = regardless of A 
or B, implying 

P{x,y)^h{x + iv) (23) 

where h is some arbitrary function. 

Consider first m = 0. Equations (fT9|) and ((2T|) then have two possible 
solutions 

A = (mdW'{z)^W'_^{z) = 



B = a.ndW'{z)^W'_{z)^ 



-E 
v//(^).9W 



corresponding to motion away from (+) and toward (-) the horizon, chosen to 
be at z = 1/a. 

Since the solution [A, 0,0,0] is an eigenvector of 7^ and has a negative 
eigenvalue its spin and momentum vectors are opposite, which is consistent 
with the fact that the particle is moving toward the horizon and the spin is up. 
The solution [0, 0, B, 0] is also an eigenvector of 7^ with positive eigenvalue; its 
spin and momentum vectors are therefore in the same direction, consistent with 
the particle being spin up and moving away from the horizon. 

Hence with the Rindler horizon at z = 1/a the (±) cases correspond to 
outgoing/incoming solutions of the same spin. Note that neither of these cases 
is an antiparticlc solution since we assumed positive frequency modes as a part 
of the ansatz. In computing the imaginary part of the action we note that 
P{x,y) must be complex (other than the trivial solution of P = 0), and so will 
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yield a contribution. However it is the same for both incoming and outgoing 
solutions, and so will cancel out in computing the emission probability 



Pi-oh[out] cxp[-2(Im W+ + Im h)] 



(24) 



Prob[Mi] exp[-2(Im W_ + Im h)] 
exp[-2(ImM^+ -lmW^)= exp[-4ImVF+] (25) 



using reasoning similar to the scalar case. We obtain 

W+{z) 



and after integrating around the pole (and dropping the + subscript) 

niE iriE 
W = , = — 26) 

The resulting tunnelling probability is 

r 27r , 
r = exp[ — —E] 

recovering 

^ ^ (27) 

which is the Unruh temperature. 

In the massive case equations (fT9|) and ([21]) no longer decouple. We will 
start by eliminating the function W'{z) from the two equations so we can find 
an equation relating A and B in terms of known values. Subtracting B x pip 
from Ax (fT9)l gives 

2 ABE ,^ ^ 

, + mA^ - mB^ = 



m^/m{-)^ + 2E{-)-rny/m = 



and so 



where 



A _ -E± y/E^+m^f{z) 



lim M±VE^m] . 



z^zo 



lim 

z^zo 



-E- y^E^ +m^f{z) 
my/J{zj 
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Consequently at the Rindler horizon either ^ ^ or ^ ^ —oo, i.e. either 



A ^ or i? — > 0. For A ^ at the horizon, we solve ((2T|) in terms of m and 
insert into (|T9l) 

-E ,, \ ( E 



W'{z) = W^l(z) ^ 



whereas for B — > at the horizon we solve p9p in terms of m and insert into 



2Tj) to get 



W'{z) = W'_{z) = 



-E ll + f^ 



Since the extra contributions vanish at the horizon, the result of integrating 
around the pole for W in the massive case is the same as the massless case and 
we recover the Unruh temperature for the fermionic Rindler vacuum. 

The spin-down case proceeds in a manner fully analogous to the spin-up 
case discussed above. Other than some changes of sign the equations are of the 
same form as the spin up case. For both the massive and massless cases the 
Unruh temperature (j27[) is obtained, implying that both spin up and spin down 
particles are emitted at the same rate. 



4 Black Hole Fermion Emission 

We turn next to a general static spherically symmetric black hole. As stated 
in the introduction, we will ignore any change in the angular momentum of the 
black hole due to the spin of the emitted particle. This is a good approximation 
for black holes of sufficient mass. The zero angular momentum state is main- 
tained because statistically as many particles with spin in one direction will be 
emitted as particles with spin in the opposite direction. 

We will now extend the fermion tunnelling approach to a general black hole 
with spherical symmetry. The metric is 

= -f{r)dt'^ + "TT + '^'^ide'^ + sin2(6l)d(/)^) (28) 
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where for this case we wih pick for the 7 matrices 

1 



7 

1' 



i 

-i 







1 / 



a' 



r sm ( 



a2 




where we measure spin in terms of the r-direction. The matrix for 7^ is 



7^ = ij^Yl^l'^ = H 



I gir) 1 

/(r) sin ( 



-1 

1 



The spin up (i.e. +ve r-direction) and spin down (i.e. -ve r-direction) 
sohitions have the form 



V'|(t,r,0,(^) = 



A(t,r,0,</.)e^ 

A{t,r,e,c^) ' 




C(i,r,0,0)C^ 


C{t,r,e,cj)) 




exp 



exp 



exp 



exp 



:Il{t,r, 9,. 



(29) 



(30) 



We will only solve the spin up case explicitly since the spin-down case is 
fully analogous. Employing the ansatz (j29p into the Dirac equation results in 



iA 



iB 



dtl^ + By^drl^ ] +Am = 



B 

r 



del] + ^«a^/T = 



sm 



dtl] - Ay/^drl] ] +Bm = 



A 
r 



deh 



smf 



(31) 
(32) 
(33) 
(34) 
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to leading order in h. We assume the action takes the form 

/t = -Et + W{r) + J{e, (j>) (35) 

where we only concern ourselves with positive frequency contributions as before. 
This yields 

/ iAE 



.VKr 



iBE 



B^/g{r)W'{r)\ +mA = (36) 



-(jo + -^iJc^ ] = (37) 



Ay/^W' (r) ] + Bm = (38) 



--{Je + ^^J^ ] = (39) 



Notice that ((37|) and (|39p result in the same equation regardless of A or _B 
(i.e. ( Je + -^^iJrti) — must be true), implying that J(0, (j)) must be a complex 
function. As with the Rindlcr case, the same solution for J is obtained for both 
the outgoing and incoming cases. Consequently the contribution from J cancels 
out upon dividing the outgoing probability by the incoming probability as in eq. 
([7]). We therefore can ignore J from this point (or else pick the trivial J = 
solution) . 

Equations and (|38p (for to = 0) have two possible solutions: 
A ^ -iS and VF'(r) = W+(r) ^ 



A = iB a.TLdiW'{r) =W'_{r) ^ 



~E 



where W+ corresponds to outward solutions and W- correspond to the incoming 
solutions. The overall tunnelling probability is 

Prob[o^t] exph2(ImTy+)] 

i = — n — r = f — ttt^ — ,, — exp — 4im W+\ (40) 

Prob[m] exp[-2(ImVF_)] ^ ^ ^ ' 

with 



After integrating around the pole (and dropping the + subscript) we find 

W = ^ (41) 
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giving 

for the resultant tunnelling probability to leading order in h. 
We therefore recover the expected Hawking Temperature 



Th ^ (43) 

in the massless case. 

Solving equations (f36|) and ([38|) for A and B in the case that m ^ leads to 
the result: 

■aV _ ~iE+^f(^m 



and approaching the horizon we see that linir^ro (5)^ = ~1- Following a 
procedure similar to what was done above, we obtain the same result for the 
Hawking Temperature as in the massless case. 

The spin-down calculation is very similar to the spin-up case discussed above. 
Other than some changes of sign, the equations are of the same form as the 
spin up case. For both the massive and massless spin down cases the Hawking 
temperature (I43p is obtained, implying that both spin up and spin down particles 
are emitted at the same rate. This is consistent with our initial assumption 
that there are as many spin up as spin down fcrmions emitted. 

4.1 Painleve Coordinates 

In this section we demonstrate that Painleve coordinates can be used to recover 
the results of the preceding section, albeit by a somewhat different computa- 
tional route. 

Using the transformation 




we obtain from the metric 

ds^ = -f{r)dt^ + 2y/J(r)J-^ - Idrdt + dr^ + r^dfl^ (45) 



which is the Painleve form of a spherically symmetric metric. 

This coordinate system has a number of interesting features. At any fixed 
time the spatial geometry is flat. At any fixed radius the boundary geometry 
for the Painleve metric is exactly the same as that of the unaltered black hole 
metric. Also, this form of the Painleve metric is a very convenient metric to use 
for black hole tunnelling since the imaginary part of the action for the incoming 
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solution is zero which means Prob[m] — 1 [T^. This property also holds for 
fermion tunnelling. 

We choose the representation for the 7 matrices to be 



7* 



7 



^ = r 



l + ^l-g{r)a 



rr^f a3 



1/0 ai 



1 / o- 



7 = — ^ ^2 



r sm ( 



a' 



The matrix for 7^ for this case is: 



7^ = ij^Yl^l^ 



1 



/(r) sin ^^ 



-1 - ^/i - .gWc^^ 

+1 - Vl - 9ir)<J^ 



Measuring spin in the r-direction we have as before the two following ansatz for 
the spin 1/2 Dirac field which correspond to the spin up (i.e. +ve r-direction) 
and spin down (i.e. -ve r-direction) cases: 



V'|(t,r,0,(^) = 



A(t,r,0,</))CT 

A{t,r,e,c^) 






C{t,r,B,4>) 


D{t,r,0,(t)) 



exp 



exp 



exp 



exp 



(46) 



(47) 



Once again we will only solve the spin up case explicitly. Insertion of the 
ansatz into the Dirac equation results in the following equations to the leading 
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order in h. 



B ( ^= (l + - 9ir)) dtl^ + VsM^./t I + ^"^ = (48) 



- (del^ + -^td^I^ ] = (49) 



J= (l - Vl-5W) dtl^ - TRO^r/T ) + = (50) 



- ( dgiT + -^id^lT ] = (51) 
r \ smO ' 



x/7W 



To solve these equations we pick the ansatz ([35]) for the action, again working 
only with positive frequency contributions. The equations for J are the same as 
in the last section, and we can dispense with this function for the same reasons 
as before. We obtain 



B 



= (l + - 9{r)) E - ^W)W'{r) J + Am - (52) 



= (l - - 9ir)) E + ,/WW'{r) \ + Bm = (53) 



Equations and (|53p (for to = 0) have two possible solutions: 



E 

A = and W'{r) = W'_^{r) = — 



B = Q&uAW'{r)=WL{r) = 



(l + Vl-5(0) 



-El 



W+ corresponds to outward solutions and W- correspond to the incoming so- 
lutions. Notice that have a pole at the horizon but W'_ has a well de- 
fined limit at the horizon and does not have a pole (i.e. lim^-^ro W- {'"') = 

\J /'(''o) ■ ^^^^ implies that the the imaginary part W_ is zero and confirms 
that Prob[m] — 1. So the overall tunnelling probability is: 

r oc Prob[oMt] 

r (X exp[-2ImVK+] 



W+{r) = 



e{i + ^/T^W)) 



dr 
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and after integrating around the pole (and dropping the + subscript): 

W = ^ (54) 
V5'(r-o)/'(ro) 

So the resulting tunnelling probability is once again: 

47r 

r = cxp[ E] 

and the normal Hawking Temperature is also recovered for the Painleve massless 
case 

Th = ^ (55) 

Solving equations and ([55)1 for A and B in the case that m ^ leads to 
the results that A ^ as r — > rp or _B ^ as r — > tq . So the same final result 
will be recovered in the massive case. 

4.2 Kruskal-Szekers Metric 

In the preceding subsections we employed metrics that had co-ordinate singu- 
larities at the horizon. One might be concerned that the tunnelling effect 
is dependent upon this. Here we demonstrate that this is not the case, by 
investigating fermion tunnelling in the Kruskal-Szekers metric 

ds^ = f{r) {-dT^ + dX^) + r'^dn^ (56) 

where: 

■'^ ' r 4m ' 



The metric ([56| is well behaved at both the future and past horizons X — ±T 
(corresponding to r = 2M). Note that the metric has a timelikc Killing vector 
Xdr + Tdx (and not dr)- 

For this calculation we will employ the following representation for the 7 
matrices 



7^ = 




' 


1 









7^ = 




' 











1' = 


1/ 
T \ ai 





) 


1" = 


' { 







r sin \ 


a2 






where we measure spin referenced to the X-direction. The matrix for 7 is 
75 = z7*7'-7«7^ ^ 1 / -1 



/(r) r2 sin 6* V 1 



15 



The spin up (i.e. +ve X-direction) and spin down (i.e. -ve X-dircction) 
solutions have the form 



■<Pl(T,X,y,z) 



A{T,X,e,cj))^^ 
B{T,X,e,4>% 
A{T,X,e,cp) ' 


B{T,X,e,cp) 


c{T,x,e,4,)i^ 
D{T,x,e,4>)i^ 



c{T,x,e,cp) 



D{T,X,e,cp) 



exp 



exp 



exp 



exp 



^/t(t,x. 



j:Il{T,X,e, 



-IliT,X,9,. 



(57) 



(58) 



Once again inserting the spin-up ansatz (|57p (the spin-down case being similar) 
into the Dirac equation yields the following equations 



B 



[drl] +dxI]) + A'm 



1 



A 



A 



{dxl] ~ dxl] ) + Brn 
1 



del 



sm( 



(59) 

(60) 
(61) 

(62) 



to leading order in h. This time we can infer only that the action takes the form 

/t = -/(x,r) + J(^,<^) (63) 

The equations for J are unchanged from previous calculations. We thus ig- 
nore these equations since they do not affect the final result and only concern 
ourselves with solving for I[X,T). 

In order to solve the equations we need a definition of the energy of the wave. 
We will define energy via the timelikc killing vector 

= NiXdr + Tdx) 

where TV is a normalization constant chosen so that the norm of the Killing 
vector is equal to 1 at infinity. This yields 

dx^^iXdr + Tdx) (64) 

and so 
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d^I = -E (65) 

Using ([55)1 with ([55]) and (pT|) we shall solve the equations. 

Consider first the massless case. Here either A = Q or B — Q. For A — 
(outgoing case): 

drl + dxl = 
±.^XdTl + TdxI) = -E 

The first equation implies the general solution of / = h{X — T) and the second 
in turn leads to 



AME = {X - T)h'{X - T) 
AME 



h'{X - T) 



[X-T) 



which has a simple pole at the black hole horizon X = T. Setting rj = X — T 
we have 

V 

Integrating around the pole at the horizon (doing a half circle contour) 
implies 

ImIout=ATrME 

for outgoing particles. 

For the incoming case B — and so 



drl - dxl = 
-l-^{XdTl + TdxI) = -E 

The first equation implies the general solution / = k{X + T) and so the second 
leads to 



-AME = {X + T)k'{X + T) 
-AME 

+ - JxTt) 

Note that this equation does not have a pole at the black hole horizon X = T. 
Hence for incoming particles 

Im lin = 
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and so Prob[m] = 1 like in the Painleve case. The final result for the tunnelling 
probability is 



Prob[out] 



exp[-2Im/o„t] = exp[-87rAf£'] 



Prob[m] 

and we see that the Hawking Temperature Th ~ s^m recovered in the mass- 
less case. 

In the massive case we must use equations ((65|). ([59l) and (|6T1) to solve for 
^ . A straightforward calculation yields 



(67) 



where we note as the black hole horizon {X = T) is approached that either 
I ^ or I ^ = Subtracting ^/A from ^/B 



leads to 



y/f{2M)vaT y/f{2M)mX ' 



drl = —dxl 



(i + (4)^) 



and so from (j65p we obtain 
dxl = 



[x(i-(|)2) + r(i + (|)2)] 



From eq ((67)) we find that 



where ^ ^ at X 



(68) 



lim 



and 

^'^TdX 



= lim 



.A, 



A^ d ,A, 



(l-(^)^) + 2(X + T)(-)-(-) 



= 1 



Consequently dxl has a simple pole at the black hole horizon implying Im lout = 
AttME in the massive case. Not that when 4 , ^^^'^ then dxl does 



B ^f{2M)mT 

not have a pole at the horizon, implying that Im — 0. The rest of the 
calculation proceeds as before, and we recover the Hawking temperature in the 
massive case. 



4.3 Conclusions 

We have shown for the first time that computing the Unruh and Hawking tem- 
peratures using the tunnelling method holds for fermions. Comparatively few 
demonstrations that fermionic radiation has the same temperature as scalar 
radiation due to the presence of these horizons appears in the literature [24] - 
[27] . These all involve either lower dimensional calculations of the Bogoliubov 
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transformation [21] or use of the GTCT [IS]- [57] to calculate fermion radiatfon 
from evaporating black holes. We have shown that computation of black hole 
temperature for fermion emission using tunnelling methods is relatively simple 
and straightforward to compute. 

For accelerated observers using Rindler coordinates we found the expected 
Unruh temperature. We also applied fermion tunnelling to a general static 
spherically symmetric black hole metric in both Schwarzschild and Painleve 
form, and found that the usual Hawking temperature is recovered. That 
this situation does not depend on coordinate singularities was demonstrated 
by showing the same results hold for the Kruskal-Szekers metric. Our results 
indicate not only that the tunnelling method is robust, but that it can indeed 
be understood as a physical phenomenon. 

Extending fermion tunnelling to rotating spacetimes in which the emitted 
particles have orbital angular momentum would be a natural next step. Com- 
puting fermion tunnelling in the background of the Kinnersley metric is a nat- 
ural step. Based on the emission probability results from the Kinnersley Black 
hole [57] , we expect that the final tunnelling probability should be of the form 
exp{—^{E—V,HJ^ + C)), where C parametrizes the coupling between the spin 
of the field and the angular momentum of the black hole. Extending fermion 
tunnelling to dynamical black holes such as Vaidya or those used in [T7] would 
be a logical next step. Computing corrections to the tunnelling probability by 
fully taking into account conservation of energy will yield corrections to the 
fermion emission temperature. In various scalar field cases this is inherent in 
the Parikh/Wilczek tunnelhng method [3], [I]-[T8] and can be incorporated into 
the Hamilton- Jacobi tunnelling approach [6]. Another avenue of research is 
to perform tunnelling calculations to higher order in WKB (in both the scalar 
field and fermionic cases) in order to calculate grey body effects. The possibility 
of calculating a density matrix for the emitted particles via the tunnelling ap- 
proach in order to calculate correlations between particles is another interesting 
line of research. Work on these areas is in progress. 
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